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Abstract 

In this article we propose a dynamic approach to complex vector reconstruction in the con¬ 
text of quantum tomography. There are two underlying assumptions behind our reasoning. The 
first one claims that the evolution of a d-level pure quantum system is given by the Schrodinger 
equation with a time-independent Hamiltonian and the other states that the knowledge about 
the quantum state is provided from projective measurements, also called intensity measurements. 
The problem of quantum state reconstruction is connected with the notion known as phase re¬ 
trieval - recovering a complex vector from modulus of inner product with frame vectors. Phase 
retrieval is widely studied in many areas of science but still there is a number of problems that 
remain to be answered. We believe that the dynamic approach can significantly improve the 
effectiveness of the vector reconstruction as it aims to decrease the number of distinct projectors 
by taking advantage of the knowledge about the evolution. General conditions and observations 
are applied to a specific unitary evolution model. 

PACS: 03.65.Wj, 03.65.Aa, 03.65.Ta 

Keywords: quantum tomography, phase retrieval, state reconstruction, complex vector recon¬ 
struction, intensity measurements 


1 Introduction 

To clarify the problem, let us postulate that the achievable information about a d-level pure quantum 
system is encoded in a complex vector - called the state vector and denoted which belongs to 
the Hilbert space H = C d such that dirriH = d < oo. Moreover, \ip) is normalized, i.e. (i/j\ ip) = 1, 
where (.|.) denotes the inner product in the Hilbert space. Then the goal of quantum tomography is 
to reconstruct the accurate representation of \ip) on the basis of data accessible from an experiment. 
Naturally, multiplying the state vector by a scalar of unit modulus does not change the measurement 
results. Thus, the state vector can be determined up to a global phase factor. 

Apparently, there are many approaches to quantum tomography which differ from one another 
when it comes to the kind of measurement(s) and the number of their repetitions. One can refer 
to papers [1, 2, 3, 4, 5] to get some insight about the methods used in quantum tomography. In 
case of pure states tomography, we are analyzing the problem of recovering a complex vector from 
intensity measurements - the very same kind of problem is considered in many other areas of 
science, from pure mathematics to speech recognition or signal processing. Thus, there is a vast 
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literature concerning phase retrieval. In recent years the attention was paid by many researchers 
to the connection between complex vector reconstruction and the theory of frames - out of many 
papers one can especially refer to [6, 7, 8, 9, 10, 11]. By an N —element complex frame in C d , 
denoted <f> = ..., |0 jv)} (where 1 9i) € C d ), one should understand a set of complex vectors 

that span C d . In articles not connected to quantum tomography authors usually consider in general 
the problem of reconstructing an unknown complex vector x € C d from its intensity measurements, 
i.e. it is discussed whether the knowledge about the non-linear map given by 

Jq> : X —>• (|(#i|£)| 2 )i=l,...,JV) (1) 

is sufficient to determine the complex vector x. However, in papers not related to Physics authors 
do not usually discuss the problem of how the intensities can be obtained. Whereas in quantum 
tomography there is a tendency to look at the reconstruction problem from the point of view of 
’economy of measurements’. Since each distinct kind of measurement requires, in general, preparing 
a new experimental setup, quantum physicists would prefer to perform quantum tomography by 
the lowest possible number of distinct kinds of measurements. It turns out that if one knows 
how the quantum state changes in time, one can obtain new data by performing the same kind of 
measurement at different time instants (apparently, each time one has to measure a distinct physical 
system). This observation gives the gist of the stroboscopic (dynamic) approach to complex vector 
reconstruction. 

We believe that the stroboscopic approach seems the most advantageous as it focuses on de¬ 
termining the optimal criteria for quantum tomography. This approach originated in 1983 in the 
article [12], in which the author was considering the minimal number of distinct observables for 
quantum tomography of a system with evolution given by the von Neumann equation. Later, the 
stroboscopic tomography was developed in many papers, from which one can especially refer to 
[13, 14, 15]. Recently new results concerning the stroboscopic tomography have been proposed in 
[16]- 

Originally, the stroboscopic tomography was applied to mixed quantum states with evolution 
given by the Kossakowski-Lindblad equation of the form pt = L[pt], where the operator L is called 
the generator of evolution [17, 18]. In order to reconstruct the initial density matrix po one has 
to assume that it is possible to obtain from an experiment mean values of certain observables 
from a fixed set {Qi,. .., Q r }, which does not satisfy the condition for completeness. Provided one 
can measure the same observable more than once there may exist a set of discrete time instants 
{ti,..., ip}, such that the matrix of data (Qi)tj = Tr(Qip tj ) where i = 1,..., r and j = 1 ,p 
accessible from an experiment enables one to reconstruct the initial density matrix. Moreover, 
knowledge about the evolution makes it possible to determine not only the initial density matrix, but 
also the complete trajectory of the quantum system. The minimal number of distinct observables 
needed for quantum tomography as well as the algebraic structure of those observables can be 
determined on the basis of the generator of evolution. Explicit theorems concerning density matrix 
reconstruction in this approach can be found in [15]. 

In this article we reformulate the assumptions of the stroboscopic tomography so that it will be 
applicable only for pure states with unitary evolution. One might think that it will be a simplification 
of the original stroboscopic approach as introduced in [12]. However, one should bear in mind that 
in case of density matrix reconstruction there are specific general criteria for quantum tomography 
(concerning the minimal number of observables, time instants etc.), whereas in case of complex 
vector reconstruction there are still many open questions - for instance it remains unknown what is 
the minimal number of intensity measurements for phase retrieval (see more in section 2). Therefore, 
this paper is not just an adaptation of the stroboscopic tomography to pure states, but it tackles in 
an original way current problems connected with phase retrieval, i.e. complex vector reconstruction. 
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Moreover, additional motivation behind this article is the observation that in quantum optics or 
quantum communication quite often one can be sure that the quantum state of a physical system 
is pure and, therefore, the achievable information is encoded in a complex vector. Thus, we believe 
that in such situations it is more efficient to employ a quantum tomography approach devised 
exclusively for pure states. Furthermore, we believe that the reasoning introduced here might be 
applicable to phase retrieval problems that arise in other areas of science provided one can access 
knowledge concerning how the vector changes in time. 

In this article we consider the case when the evolution of a d— level quantum system is given by 
the Schrodinger equation of the form 

- mm), ( 2 ) 

where | ip(t)) € Ti and H is a self-adjoint operator such that H : H —>• H. The symbol h denotes the 
Dirac’s constant, for simplicity, henceforth, it is assumed that h = 1. Moreover, the vector space of 
all linear operators on % shall be denoted by 13(H). 

If the evolution of the state vector is given by (2), then \ip(t)) for any t € M + can be computed 
from 

\m) = umm, m 

where U(t) € 13(H) is a unitary operator such that U(t) = exp(-iHt). 

Discussing the problem of quantum tomography, we will assume that the knowledge about the 
quantum system is provided by results of projective measurements. Thus, we assume to have a 
set of operators {Mi,..., M r } for some r G N, and each Mi € 13(H). The operators are assumed 

to be projectors, i.e. Mi = |i)(z|, where |i) € H for * = 1,..., r. Naturally, the vectors |z) are 

normalized, i.e. (i|i) = 1 for i = 1,... ,r. We assume that the set of projectors is incomplete and 
single measurement of each projector does not allow one to reconstruct the initial state vector |-0(O)). 
Thus, in order to determine the initial state vector one has to perform each projective measurement 
at some discrete time instants {t±,... ,tp}- Consequently, the data obtainable from an experiment 
can be expressed as rrii(tj) = ('i/j(tj)\Mi\'ip(tj)). In section 2 we consider how such data can be 
used for effective complex vector reconstruction. Obviously, we assume to have a large number of 
identically prepared quantum systems and, therefore, each individual system is measured only once, 
which allows us to skip the problem of how the measurement changes the state. 

In this section we have formulated the problem of quantum tomography for pure states and enu¬ 
merated the assumptions of the stroboscopic approach. In section 2 we shall propose general results 
and observations that are obtainable from employing the assumptions and combining them with 
current knowledge about stroboscopic tomography and phase retrieval. In section 3 we demonstrate 
how the stroboscopic tomography works for pure states by solving a specific example. 

2 General Results of the Stroboscopic Tomography for Complex 
Vector 

First, we shall expand the exponential formula for the unitary operator U(t) introduced in (3). The 
most obvious way to expand this operator goes as follows 

U(t) = exp(-iHt) = J2 [ - Z y L n k - ( 4 ) 

k =0 

However, one can agree that it should be possible to expand exp(—iHt ) by means of a finite sum 
of powers of H due to the fact that for every matrix one can introduce the notion of the minimal 
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polynomial. If by p one denotes the degree of the minimal polynomial of H, then the decomposition 
(4) can be rewritten in the form 


M-i 

U(t) = exp(-iHt) = oik{t)H k , (5) 

k =o 

where atk(t) are certain time-dependent functions. It can be proved that they are mutually lin¬ 
early independent and can be computed from a set of differential equations which depend on the 
coefficients of the minimal polynomial of H [14]. 

The exponential form of the unitary operator is useful to analyze the experimental data. The 
formula for the measurement result can be written as 


miitj) = ('i/j(t j )\Mi\'i/;(t j )) = {'ip(t j )\i){i\i>{t j )) = 
= \{i\Tp{tj))\ 2 = |(z|exp(-ifff i )|V’(0))| 2 


M-l 

k =0 

1 

k =0 


M-l 

k =0 


( 6 ) 


where 1 4^) = H k \i). 

Now, if one assumes that projective measurement associated with the operator Mi is performed 
at time instants {t \,..., t p }, we obtain the set of p equations concerning each operator Mp. 

m— i 

mi{ti) = Y M*i)| 2 |(#V(0)>| 2 , 

k =0 

M-l 

mi(t 2 ) = Y l«fe( i 2 )| 2 |(</>fV( 0 )}| 2 , 

k =o (?) 


M— 1 

m(.t P ) = Y MrfK^fVto))! 2 . 

k =0 

One can notice that such a system of equations can be rewritten as a matrix equation 




|ao(^i)| 2 l«i(ii)| 2 ••• |a^_i(ii)| 2 


' i(4°V(o))i 2 ' 

mi(t 2 ) 

= 

\a 0 (t2)\ 2 |ai(t2)| 2 ••• |a M -i(£2)| 2 


K^'W))! 2 

rrn(tp)_ 


_|a 0 (t p )| 2 \ai(t p )\ 2 ... \ap-i(t p )\ 2 _ 


JWA‘ ) l</’(o)}| 2 . 


On the left-hand side of the equation (8) we have a vector of data accessible from an experiment. 
The matrix A = [|afc(t, )| 2 ] where k = 0,..., p — 1 and j = 1,... ,p is computable on the basis of the 
minimal polynomial of H. Thus, if the matrix A is invertible, one can calculate from the equation 
(8) the intensities (projections) |(^ fc ^(0))| 2 for k = 0,..., p — 1. The condition for computability 
of these projections can be stated as a theorem. 
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Theorem 2.1. From the matrix equation (8) one can calculate the projections \(4>\ k ^\f’(Q))\ 2 f or 
k = 0,..., /a — 1 if and only if 

P = h-, (9a) 

detA / 0. (9b) 

One can easily notice that the condition (9a) ensures that the matrix A is square and the condition 
(9b) ensures its invertibility. 

Theorem 2.1 states that the number of time instants is equal to the degree of the minimal 
polynomial of H. 

Now it should be clear that if the stroboscopic approach to quantum tomography is applied to 
each projective operator from the set {Mi,..., M r }, i.e. each projector is measured at time instants 
{ti,... ,tu } such that the condition (9b) is fulfilled, then one can determine the correspondence 

IV’(O)) —»• |(</>f' ) |'0(O))| 2 for k = 0,... — 1 and i = 1,..., r. (10) 

The set of vectors that we project |^(0)) onto shall be denoted as d> = 0 )- Hence¬ 

forth, the non-linear map that assigns to the unknown state vector IV’(O)) the set of squares of 
absolute values of its inner product with the vectors from the set as demonstrated in (10) shall 
be denoted by J$>. 

The natural question which arises here states: When can we reconstruct the initial state vector 
on the basis of intensities defined by the set $ ? 

Apparently, of one wants to reconstruct |^(0)) on the basis of the intensities given by the 
correspondence (10), the vectors that belong to have to span the Hilbert space to which |^>(0)) 
belongs. In other words, the necessary condition for complex vector reconstructability claims that 
the collection of vectors has to constitute a frame in H such that IV'(O)) € H. 

In order to formulate the theorem concerning the necessary (but not sufficient) condition for 
reconstructability of |^(0)) let us first introduce the notion of the Krylov subspace. 

Definition 2.1. The Krylov subspace, which shall be denoted by K^(H, Mf), is defined as follows 

Kfi(H, Mi) := Span{\^), |^ 1} ),..., |^ _1) >}. (11) 

The stroboscopic approach to tomography is the more advantageous, the higher the dimension 
of the Krylov subspace is. One can instantly notice that the dimension of the Krylov subspace 
depends on algebraic properties of both |i) and H. 

The theorem concerning the necessary condition for phase retrieval can be formally formulated 
as follows. 

Theorem 2.2. The necessary condition for the initial state vector |?/>(0)) G H to be reconstructible 
on the basis of the intensity measurements |((^ fc ^|'0(O))| 2 for k = 0,...,//— 1 and i = 1,... ,r claims 
that 

r 

Q)K fl (H,M i )=H, (12) 

i =1 

where he symbol 0 denotes the Minkowski sum of subspaces. 

To consider the sufficient condition for complex vector reconstruction let us revise the general 
considerations concerning the question when phase retrieval is possible. In a recent paper [10] the 
authors claim that phase retrieval is possible when any two vectors |?/>(0)) and |'<//(())} with identical 
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intensity measurements differ only by a scalar of norm one, i.e. |^(0)) = e ie \ij/(0)). In other words, 
the same observation can be stated that it is possible to reconstruct the complex vector (0)} if 
and only if the non-linear map is injective and is a frame. Thus, henceforth in situations 
when phase retrieval is possible, we shall say that the frame $ generates (or defines) injective 
measurements. 

In [9] Bandeiraa, Cahill, Mixton and Nelson postulated a conjecture according to which if one 
wants to reconstruct a vector x such that x € C d , then a frame that contains less than Ad — 4 
vectors cannot generate injective intensity measurements, i.e. according to the authors fewer than 
Ad — 4 modulus of inner product of x with other vectors is not sufficient to obtain the structure of 
x. Furthermore, in the same paper the authors postulated the second part of the conjecture that a 
generic frame with Ad — A vectors (or more) generates injective measurements on C d . The second 
part of the conjecture has recently been proved in [10], where the authors explained the notion of a 
generic frame and proved that for a generic frame d> that contains at least Ad —A elements the map 

is injective. 

Another recent paper [11] proves a result that contradicts the first part of the conjecture from 
[9]. Cynthia Vinzant proposes a frame in C 4 which consists of 11 vectors and proves that it defines 
injective measurements on C 4 . Therefore, the current knowledge about the phase retrieval problem 
does not give an answer to the question what is the minimal number of elements of the frame <h so 
that the map is injective, i.e. so far it remains unknown in general how many intensity measure¬ 
ments at least are needed to reconstruct the unknown complex vector. However, in [9] the authors 
propose a relatively efficient way to check whether a frame generates injective measurements. 
Their approach is presented below as a theorem. 

Theorem 2.3 (Bandeiraa et. al. 2014). A frame = {|@i),..., |#at)} (where \Of) € C d ) defines 
injective measurements, i.e. one can reconstruct some unknown vector x £ C d from intensity 
measurements |(0j|x)| 2 for i = 1,..., N, if and only if the linear space 

U := {Q G C dxd : <0i|Q|0i> = ■ ■ ■ = (0n\Q\On) = 0} (13) 

does not contain any non-zero Hermitian matrix of the rank < 2. 

The theorem 2.3 states clearly the sufficient condition that needs to be fulfilled so that a frame 
<f> defines injective measurements and, therefore, it is possible to reconstruct a complex vector on 
the basis of its intensity measurements. For a given frame one can relatively easy check whether 
the condition stated in theorem 2.3 is fulfilled or not. However, so far there have been no proposal 
concerning the procedure how to obtain such a sufficient frame. 

3 An Example of the Stroboscopic Tomography Model for a Com¬ 
plex Vector from % = C 2 

In this section we demonstrate how to apply the general results and observations concerning complex 
vector reconstruction to a specific quantum tomography problem. 

Let us analyze a quantum system associated with the Hilbert space H such that dimH = 2. 
The knowledge about the system is encoded in a time-dependent state vector | ip(t)) G H. The 
problem of reconstructing the initial state vector |^(0)) from experimentally accessible data has 
been widely studied. In [19] one can find the result according to which the initial state vector can 
be uniquely determined from projections onto 4 complex vectors which correspond to the vertices of 
a tetrahedron embedded within the Bloch sphere. The main goal of this section is to demonstrate 
an alternative approach to qubit reconstruction - by employing the stroboscopic tomography one 
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can determine |^(0)) on the basis of two distinct projectors, but each o them is measured at two 
different time instants. 

As the stroboscopic approach to quantum tomography requires the knowledge about the evolu¬ 
tion of the system, let us postulate that the state vector evolves in accordance with the equation 

& = -mm), (Hi 

where the Hamiltonian, in this example, takes the form 

H = '» = (? o’) ■ < 15 > 

The goal of the quantum tomography is to determine the initial state vector |'i/>(0)) form certain 
intensity measurements. It is commonly known that the structure of a two-level state vector can 
be represented by two real parameters 6 and (j). Thus, in order to reconstruct |^(0)) one needs to 
fin 6, <f> such that 

i*°» - U&) ■ < i6) 

where 9 € [0,7r] and (j) € [0, 2tt\ [20]. 

One can notice that for H given by (15) the degree of the minimal polynomial is equal to the 
degree of the characteristic polynomial - both are equal 2. Bearing in mind representation (3) on 
can write 

M*)> = exp{-iHt) |^(0)), (17) 

where it is convenient to substitute H = — iH , which on the basis of (5) gives us 

I i'it)) = exp(Ht)\il)(0)) = (a 0 (f)Il 2 + ai(t)#) IV’(O)). (18) 

The functions ao(t) and a±(t) can be computed from the set of differential equations that they 
need to fulfill [14]. One can easily get 

ao(t) = cos(t ) and aq(t) = sin(t). (19) 

Thus, |V’(£)) can be expressed as 

\ip(t)) = cos(t)|V’(0)) + sin{t)H\i’( 0)). (20) 

To determine the initial state vector |^(0)) we propose to perform the projective measurements 
onto two vectors |1) and |2) given by 

|1J = 7S ("^ “ d |2> = 775 (?) - (21) 

In other words, we introduce projectors M\ = |1)(1| and M 2 = |2)(2|. Each distinct projective 
measurement is performed twice, as the degree of the minimal polynomial is equal 2. Let us denote 
the time instants by t\ and £ 2 - 

One can easily calculate H\l) and H |2) and obtain 

d|1> = 71 (-1) and ^ |2 > = 75 (2)' (22) 
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Obviously, the four vectors introduced in (21)-(22) constitute a frame in C 2 . 

As each projective measurement is performed two times, let us denote the results of the measure¬ 
ments by rrii(tj) = \(i\^{tj))\ 2 where i = 1,2 and j = 1,2. For each projector we get two equations 
which can be combined into a matrix according to (8). For the projector Mi we obtain 


mi(ti) 


cos 2 t\ 

sin 2 ti 

\ l( 1 IV’(o))| 2 ' 

mi(t 2 ) 


cos 2 t 2 

sin 2 t 2 

U(i|^#(o))| 2 J 


For simplicity we can assume that t\ = 0. Then from (23) we can calculate the intensities if 
sin 2 t 2 i=- 0, which means that 

0 t- kn for k £ N U {0}. (24) 

Assuming that the condition (24) concerning the time instant t 2 is fulfilled, we can calculate the 
intensities 

|(l|V>(0))| 2 =mi(0), (25a) 

|(l|fl|^(0))t 2 = + ™lM (25b) 

sin z t 2 

In case of the projector M 2 one can write a matrix equation analogous to (23) which is solvable 
under the very same condition concerning the time instant t 2 . Thus, one can get the projections 


\m\m)\ z = 


I(2|^(0))| 2 =m 2 (0), 

—cos 2 t 2 m 2 (0) + m 2 (t 2 ) 


sin 2 to 


(26a) 


(26b) 


This analysis proved that by performing the two projective measurements given by M\ and 
M 2 at two different time instants 1 1 = 0 and t 2 satisfying (24), we can determine the intensity 
measurements defined by a complex frame 4> = {11), H\l), 12), H\2)}. The question that remains to 
be answered states whether this frame defines injective measurements. To answer this question we 
propose to employ theorem 2.3. Thus, we shall analyze the linear space as introduced in (13). In 
this example it takes the from 


U = {Q € C 2x2 : (1|Q|1) = (1\HQH\1) = (2|Q|2) = (2\HQH\2) = 0}. (27) 


According to the theorem 2.3 the frame 4* defines injective measurements if and only if the linear 
space £<j> does not contain any non-zero Hermitian matrix of rank <2. In the analyzed case the 
proof of the injectivity can be done by a contradiction, i.e. we shall assume that there is a Hermitian 
matrix Q' which belongs to £$. Bearing in mind dimH = 2, Q' can generally be presented in the 
form 


Q' 


fa c + id\ 
\c — id b J 


(28) 


for some a, b,c,d £ M. By employing the software Mathematica 10, one can get that Q' € £$ if and 
only ifa = 6 = c = d = 0, which means that the frame 4> defines injective measurements. Therefore, 
from the intensity measurements |(i|i7 ;i j?/>(0)}| 2 for i = 1,2 and k = 0,1 one can reconstruct the 
initial state vector |V’(0))- 

In order to observe that phase retrieval is possible let us calculate the intensities defined by 4> 
from the knowledge about the vectors (21)-(22) and the general representation of 1^(0)) (16). One 
can easily get 


|<1|V>(0))| 2 


5 — 3 cos6 — 4 cos(j)sin0 


10 


(29a) 
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\m\m)\ 2 - 5+3cose+icos4 ’ sine , 

(29b) 

(n /(o)) 2 5 + 3 cost) + 4siruj)sin6 

(29c) 

- 5 ~ + isin *‘ in0 - 

(29d) 


Comparing the equations (25a)-(26b) with the equations (29a)-(29d) one can obtain 

10 m 2 (0) - ra 2 (f 2 ) 

cosO = - — -. 

b sin z t 2 

Bearing in mind that 6 £ [0,7r], one can uniquely determine 9 by 

'10 m2(0) — m2(t2) 


8 = arc cos 


6 


sin 2 t 2 


(30) 


(31) 


To find the structure of the initial state vector IV’(O)) parametrized as (16), we also need to 
determine q i. Again by comparing the set of equations (25a)-(26b) with the equations (29a)-(29d) 
one can obtain 


sine/) 


5 rri 2 {t 2 ) — rri 2 (0) cos (2^2) — sin 2 t 2 
4 sin 2 t2sin0 


(32a) 


5 2(mi(t 2 ) - cos 2 t2mi(0 )) + m 2 (t 2 ) — rn 2 { 0) - sm 2 f 2 
4 sin6 sin 2 t 2 


(32b) 


One can agree that for cj) £ [0, 27t] the knowledge about sincj) and cos(p is sufficient to uniquely 
determine cj). Therefore, we have showed explicitly that for an isolated quantum system with 
evolution given by (14), projective measurement performed for two operators introduced in (21) in 
two distinct time instants ti = 0 and t 2 satisfying the condition (24) provides sufficient data for 
phase retrieval, i.e. the experimental data is sufficient to uniquely determine the structure of the 
initial state vector |^(0)). 


4 Summary 

The problem of complex vector reconstruction from intensity measurements appears in many areas 
of science. Thus, different aspects of phase retrieval have been deeply studied. However, there 
are still many open problems, concerning for example the optimal criteria for reconstruction. This 
article proposes in the context of quantum tomography a new approach to complex vector recon¬ 
struction, which allows to decrease the number of distinct projectors required for phase retrieval. 
The stroboscopic approach has been applied to one specific model of unitary evolution, but the 
number of problems that can be solved with this reasoning is unlimited. One can also apply the 
introduced analysis to phase retrieval problems not connected to quantum tomography provided 
one can access the knowledge concerning how the vector changes in time. 
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